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ABSTRACT 


The  use  of  differentially  deflected  fairwater  planes 

i 

to  control  submarine  rolling  is  studied.  Because  of  couDling 
between  pitch  and  roll  angles,  the  snap  roll  that  occurs  in 
a  high  speed  hard  turn  affects  the  stability  of  a  submarine 
not  only  in  the  horizontal  plane  but  also  in  the  vertical 
plane.  Direct  roll  control  was  achieved  by  making  use  of  the 
fairwater  planes  in  a  differentially  deflected  mode  such  that 
they  could  give  counter  moment  to  reduce  the  snap  roll.  A 
roll  controller  was  designed  as  a  position  and  velocity  feed¬ 
back  controller.  Controlled  roll  angle  improved  the  depth 
and  pitch  stability  of  submarine. 
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101.  Sternplane  Angle  vs.  Time.  Final  Result  With  Roll  145 

Error  Limiter.  Kl  *  3.  K2  *  10.  UCK  *  12  Knots. 

Rudder  Ordered  =*  35°.  Initial  Roll  Angle  =  0° 

102.  Sailplane  Angle  vs.  Time.  Final  Result  With  Roll  146 

Error  Limiter.  Kl  =3.  K2  =  10.  UCK  *  12  Knots. 

Rudder  Ordered  «  35°.  Initial  Roll  Angle  =  0° 

103.  Roll  vs.  Time.  Final  Result  With  Roll  Error  Limiter.  147 

Kl  =  3.  K2  =  10.  UCK  =  24  Knots.  Rudder  Ordered 
=*35  .  Initial  Roll  Angle  =  -20° 

104.  Roll  vs.  Time.  Final  Result  With  Roll  Error  Limiter.  148 

Kl  =3.  K2  =  10.  UCK  =  18  Knots.  Rudder  Ordered 
=  35°.  Initial  Roll  Angle  =  -20° 

105.  Roll  vs.  Time.  Final  Result  With  Roll  Error  Limiter.  149 

Kl  =3.  K2  =  10.  UCK  =  12  Knots.  Rudder  Ordered 
=  35°.  Initial  Roll  Angle  =*  -20° 

106.1.  Depth  vs.  Time.  Poll  Stability  Step  Test.  UCK  =*  154 

24  Knots.  Rudder  Ordered  -  35° 

.2.  Disturbance  Moment  vs.  Time 

107.1.  Pitch  vs.  Time.  Roll  Stability  Step  Test.  UCK  =  155 

24  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 

108.1.  Roll  vs.  Time.  Roll  Stability  Step  Test.  156 

UCK  =*  24  Knots.  Rudder  Ordered  »  35° 

.2.  Distrubance  Moment  vs.  Time 

109.1.  Sternplane  Angle  vs.  Time.  Roll  Stability  Step  157 

Test.  UCK  =  24  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 
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110.1.  Sailolane  Angle  vs.  Time.  Poll  Stability  Step  Test.  158 

UCK  *  24  Knots.  Rudder  Ordered  ■  35° 

.2,  Distrubance  Moment  vs.  Time 

111.1.  Depth  vs.  Time.  Poll  Stability  Step  Test.  UCK  »  159 

18  Knots.  Rudder  Ordered  *  35° 

.2.  Disturbance  Moment  vs.  Time 

112.1.  Pitch  vs.  Time.  Roll  Stability  Step  Test.  UCK  »  160 

18  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 

113.1.  Roll  vs.  Time.  Roll  Stability  Step  Test.  UCK  *  161 

18  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 

114.1.  Sternplane  Angle  vs.  Time.  Roll  Stability  Step  162 

Test.  UCK  =  18  Knots.  Rudder  Ordered  =  35° 

.2.  Distrubance  Moment  vs.  Time 

115.1.  Sailplane  Angle  vs.  Time.  Roll  Stability  Step  163 

Test.  UCK  -  18  Knots.  Rudder  Ordered  ■  35° 

.2.  Disturbance  Moment  vs.  Time. 

116.1.  Depth  vs.  Time.  Roll  Stability  Step  Test.  164 

UCK  =  12  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 

117.1.  Pitch  vs.  Time.  Roll  Stability  Step  Test.  165 

UCK  =*  12  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 

118.1.  Roll  vs.  Time.  Roll  Stability  Step  Test.  166 

UCK  ■  12  Knots.  Rudder  Ordered  =  35° 

.2.  Distrubance  Moment  vs.  Time 

119.1.  Sternplane  Angle  vs.  Time.  Roll  Stability  Steo  167 

Test.  UCK  =  12  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  Moment  vs.  Time 

120.1.  Sailplane  Angle  vs.  Time.  Roll  Stability  Step  Test.  168 

UCK  *  12  Knots.  Rudder  Ordered  =  35° 

.2.  Disturbance  vs.  Time 


13 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  his  sincere  aopreciation 
to  Professor  George  J.  Thaler  for  his  encouragement  and  ■ 
guidance  during  the  preparation  of  this  thesis. 


14 


I.  INTRODUCTION 


The  importance  of  automatic  control  of  Ditch-depth  and 
roll  of  a  turning  submarine  has  become  more  obvious  with  the 
improvements  of  high  speed  nuclear  submarines.  This  study  was 
aimed  at  smoothing  the  roll  of  a  turning  submarine  in  hicrh 
speed,  which  is  known  as  snap  roll  and  causes  the  basic  problem. 
Because  of  coupling  between  pitch  and  roll,  before  attenuating 
to  smooth  the  snap  roll,  controlling  of  pitch  and  depth  in 
high  speeds  turns  was  also  studied.  In  this  study  an  auto¬ 
matic  course  controller  was  not  considered. 

The  methods  and  techniques,  which  can  be  used  to  reduce 
and  smooth  the  snap  roll,  can  be  categorized  under  two  head¬ 
ings.  The  first  concerns  changes  in  the  naval  architectural 
characteristics  of  the  designs,  such  as  increasinq  GM  (meta- 
centric  height)  and  reducing  sail  size.  The  second  category 
involves  those  alternatives  which  make  use  of  an  automatic  ship 
control  system.  Rudder  sequencing  and  speed  reduction  fall 
under  this  classification.  Reference  1  investigates  the  affect 
of  increasing  GM,  reducing  sail  size  and  speed.  In  Reference  2, 
Stamps  designed  on  automatic  roll  controller  which  makes  use  of 
rudder  sequencing  as  a  function  of  aooroach  speed  and  instan¬ 
taneous  roll  angles. 

In  recent  studies,  which  were  made  by  Naval  Ship  Research 
and  Development  Center  (NSRDC) ,  it  was  proposed  to  enhance  the 
control  of  roll  in  high  soeed  turns  by  using  the  Fairi-rater 
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planes  differentially  deflected  so  that  they  can  be  used  to 
give  a  counter  moment  to  reduce  the  roll.  The  investigation 
of  the  differentially  deflected  Fairwater  planes  effect  was 
part  of  the  project  called  "Improved  Control  For  Advanced 
Submarines".  The  project  was  carried  out  under  Program  Element 
No.  62754N  and  Task  Area  2F4 34001.  The  work  unit  number  was 
1-1563-001-74.  In  Reference  3,  estimating  the  effectiveness 
of  differentially  deflected  sailplanes  was  investigated  as  a 
part  of  project  mentioned  above.  Reference  3  was  the  inspira¬ 
tion  of  this  thesis  and  using  differentially  deflected  sail¬ 
planes  for  direct  control  of  roll  was  chosen  as  the  desion  croal 


II.  THE  NATURE  OF  THE  PROBLEM 


Turning  characteristics  of  a  surfaced  submarine,  more 
or  less,  looks  like  those  of  a  surface  ship.  But  the  situation 
in  the  submerged  position  shows  big  differences.  These  dif¬ 
ferences  are  the  result  of  different  naval  architectural  char- 
actacteristics.  Sail  structure  can  be  considered  the  main  dif¬ 
ferences  and  the  main  source  of  roll  problems.  It  is  well 
known  that  when  a  surface  ship  goes  into  a  turn,  it  exoeriences 
an  initial  inboard  roll.  After  a  very  short  transient  it  heels 
outboard.  The  reversing  of  the  roll  angle  is  primarily  due  to 
an  asymmetric  rudder.  Since  the  submarine  rudder  is  generally 
well  submerged,  the  surfaced  submarine  has  an  outboard  roll 
angle  during  both  the  initial  and  the  steady  state  phase  of  a 
turn.  For  a  typical  submarine,  at  moderate  speed  this  outboard 
roll  angle  is  less  than  10°  and  does  not  cause  any  big  problem. 

If  a  submerged  submarine  goes  into  a  high  speed  turn, 
differences  from  the  surfaced  behavior  are  noticed  and  the 
problem  becomes  three-dimensionalized.  In  the  very  first  phase 
of  the  turn,  the  ship  has  a  small  lateral  velocity  v  and  a 
small  rotational  velocity  r.  For  a  starboard  turn,  at  some 
point  along  the  centerline,  XI,  the  lateral  velocity  due  to 
r,  rXl,  equals  the  lateral  velocity  -  v.  (For  symbols  see 
Figure  1,  2,  and  Appendix  A).  The  point,  at  which  these  two 
velocities  are  equal  to  each  other,  is  called  the  instant 
center  of  turn.  At  some  point,  forward  of  the  instant  center, 
*2 ,  the  lateral  velocity,  -v,  gives  a  velocity  component  from 
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the  starboard  side.  Along  with  the  ahead  velocity  component 
of  the  ship,  these  velocities  contribute  an  angle  of  attack 
from  starboard.  At  this  point,  if  the  submarine  has  an  append¬ 
age  (sail) ,  a  lift  force  is  produced  on  the  sail  directed  to 
port.  Depending  on  the  sail  configuration  and  the  large  moment 
arm  of  the  sail.  This  force,  directed  to  the  port  gives  a  small 
outboard  (port)  roll.  In  many  submarines  this  outboard  roll  can 
not  be  felt. 

In  the  steady  state  phase  the  angle  of  attack  shifts  to 
the  other  side  (port)  and  results  in  an  inboard  roll.  But  con¬ 
trary  to  the  outboard  roll  this  inboard  roll  is  quite  signifi¬ 
cant  and  at  high  speeds  it  can  exceed  30°  which  is  considered 
hazardous  to  both  men  and  equipments.  This  inward  roll  is 
called  the  SNAP  roll.  In  Figure  .  5,  the  roll  characteristics 
of  the  ship  to  a  constant  35°  left  rudder  angle  at  24  knots  is 
shown.  As  can  be  seen,  the  initial  outboard  roll  (in  the 
simulation,  outboard  is  the  starboard  side)  is  very  small, 
less  than  1°.  But  the  snao  roll  reaches  37°. 

Snap  roll  is  not  the  only  problem  that  can  be  faced  when 
the  submarine  goes  into  a  turn.  Because  of  the  appendage 
(which  is  mainly  the  sail) ,  in  the  submerged  position,  sub¬ 
marines  do  not  have  very  good  hydrodynamics  from  this  body 
structure.  When  the  ship  experiences  a  submerged  turn  the 
sail  mainly  destroys  the  waterflow  around  the  hull.  As  a 
result  a  pressure  difference  is  created  between  the  upoer 
portion  of  the  hull  and  the  bottom  portion  of  the  hull. 
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And  a  normal  force  results  directed  downward  which  bodily 
pushes  the  ship  down.  It  is  recalled  that,  in  a  turn,  after 
the  initial  phase,  the  submarine  has  an  inward  roll  with  the 
deflected  rudder.  Rolling  inward,  causes  this  deflected 
rudder  to  act  like  a  sternplane  which  gives  the  submarine  a 
down  pitch  movement.  All  of  these  effects  come  together  and 
cause  the  ship  to  rise  or  dive  depending  on  the  peculiarity 
of  its  design.  As  it  is  stated  in  Reference  4,  "The  exact 
mechanism  of  this  is  not  well  understood".  In  Figures  3,  4, 
and  6,  the  depth,  oitch,  and  rudder  response  of  the  ship,  to 
a  35°  left  rudder  angle  at  24  knots  is  shown.  It  is  interest¬ 
ing  to  indicate  here  that  after  a  long  transient  period,  the 
submarine  starts  to  rise  with  the  bow  down  eventhough  it  dives 
initially.  It  should  be  remembered  that  in  this  simulation  of 
the  model,  no  control  surface  (stern  and  fairwater  plane)  was 
used.  The  changes  of  these  characteristics  will  be  seen  after 
the  depth-pitch  and  roll  controller  is  designed. 

The  complexity  of  controlling  the  turning  submarine 
mainly  comes  from  the  coupling  between  the  roll  and  pitch. 
Since  roll  and  oitch  are  coupled,  to  be  able  to  control  the 
roll,  which  is  our  primary  concern,  pitch  control  to  some  ex¬ 
tent  must  be  accomplished.  In  the  following  sections  a  pitch 
and  depth  controller  is  designed  by  using  only  the  sternolanes 
because  of  the  goal  of  using  the  fairwater  plane  as  a  main 
part  of  the  direct  roll  controller. 
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Figure  1.  Hydrodynamic  Forces  In  A  Turn. 
(After  Reference  4) 
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Figure  3.  Death  vs.  Time.  With  No  Controller 

UCK  *  24  Knots.  Rudder  Ordered  *  35°. 
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Figure  4.  Pitch  vs.  Time.  With  No  Controller. 

UCK  »  24  Knots.  Rudder  Ordered  *  35°. 
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X  SCALE  =  40  Seconds  Per  Inch 


Y  SCALE  =  8  Degrees  Per  Inch. 
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Figure  5.  Roll  vs.  Time.  With  No  Controller. 

UCK  »  24  Knots.  Rudder  Ordered  =*  35°. 
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Figure  6. 


Rudder  Response  vs.  Time.  With  No  Controller. 
UCK  =  24  Knots.  Rudder  Ordered  =  35°. 
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III.  DIGITAL  COMPUTER  SIMULATION  OF  SUBMARINE  MOTION 


The  mathematical  model  of  the  totally  submerged  sub¬ 
marine  is  represented  by  the  equations  of  motion  which  consist 
of  six  equations  such  that  each  of  them  represent  the  force 
and  moment  equations  on  and  around  the  three  rectangular  co¬ 
ordinate  axes.  These  equations  were  derived  by  NSRDC  in 
Reference  5  and  for  convenience  was  repeated  in  Appendix  A. 

In  this  thesis,  derivation  of  these  equations  was  not  our 
concern.  More  pronounced  knowledge  about  that  can  be  found  in 
References  4,  5,  and  6. 

The  equations  of  motion,  which  represent  the  motion  of 
the  totally  submerged  submarine,  in  six  degree  of  freedom,  con¬ 
sist  of  Hydrodynamic  coefficient,  in  undimensionalized  form 

except  the  m  (mass)  and  I  ,  I  ,  I  (three  moments  of  inertia). 

x  y  z 

In  Reference  7,  Drurey  translated  these  six  equations  into  the 
DSL  (digital  simulator  language)  form  and  originated  the  com¬ 
puter  program  which  is  the  mathematical  model  of  the  ship. 

For  the  necessity  of  using  the  hydrodynanic  coefficients  in 

undimensionalized  form,  the  equations  containing  I  ,  I  ,  I  , 

X  y  z 

both  sides  are  divided  by  l5  and  these  three  moments  of 
inertia  were  used  in  undimensionalized  form.  For  the  same 
reason  the  equations  containing  m  were  divided  by  1^. 

Since  the  trim  control  of  the  ship  was  not  considered 
in  this  study  the  ship  was  assumed  in  the  trimmed  condition. 
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For  trimming  the  ship  the  hydrodynamic  coefficient 
Z*  and  M*  were  set  to  zero.  These  two  hydrodynamic  co¬ 
efficients  represent  the  force  and  moment  acting  on  the  ship 
when  the  control  surfaces  are  at  zero  deflection. 

Since  all  of  the  simulations  of  this  study  were  done 
in  a  turn,  an  accurate  simulation  must  include  the  rudder 
actuator  dynamics  such  that  any  phase  lag  could  not  produce 
any  unstability.  Actuator  dynamics  used  in  this  thesis  was 
the  same  as  used  in  Reference  3  with  the  rudder  deflection 
angle  5°/sec.,  so  that  at  the  end  of  the  study  it  could  give 
acceptable  comparison  level  between  the  two  designs.  A  block 
diagram  of  the  rudder  actuator  is  shown  in  Figure  7. 

Since  our  prime  concern  is  to  control  the  roll  as  well 
as  the  depth  and  the  pitch  when  the  submarine  is  in  a  turn,  a 
series  of  tests  were  run  at  different  speeds  to  validate  the 
models  originated  in  Reference  7  by  Drurey.  The  model  was 
highly  satisfactory  as  shown  ir  Figures  8  through  25.  These 
figures  record  responses  in  depth,  pitch,  roll,  yaw,  speed 
changes,  and  rudder  angle  at  18,  12,  and  6  knots.  In  this 
test,  turn  was  commanded  10  seconds  after  the  simulation 
started,  to  validate  the  initial  condition  response  of  the 
model.  Since  the  submarine  was  trimmed  by  eliminating  Z*  and 
M*,  it  is  seen  that  until  the  turn  was  commanded  all  responses 
are  zero.  Simulation  at  24  knots  was  already  shown  in  Figures 
3  through  6. 

It  was  observed  that,  at  6  knots  the  ship  continues 
diving  unlike  the  simulation  of  12,  18,  and  24  knots.  The 


reason  for  this  was  that  at  low  speed  enough  pressure  dif¬ 
ference  to  cause  the  ship  to  rise  slowly  could  not  be  created. 

After  the  validation  of  the  model  it  was  decided  to  pro¬ 
ceed  with  the  designing  of  the  controllers.  In  the  following 
sections  first,  the  depth  and  pitch  controller  design  was  stud¬ 
ied.  An  optimal  control  scheme  such  as  in  Reference  7  was  used. 
The  sternplanes  was  the  control  surface  of  the  controller,  and 
the  fairwater  plane  was  saved  for  use  in  the  roll  controller. 
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Figure  8. 


Depth  vs.  Time.  With  No  Controller. 
UCK  *  18  Knots.  Rudder  Ordered  =*  35°. 
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Figure  9.  Pitch  vs.  Time.  With  No  Controller. 

UCK  =»  18  Knots.  Rudder  Ordered  *  35° 
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Figure  10.  Roll  vs.  Time.  With  No  Controller. 

UCK  *  18  Knots.  Rudder  Ordered  -  35°. 
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Figure  11. 


Speed  Change  vs.  Time.  With  No  Controller. 
UCK  ■  18  Knots.  Rudder  Ordered  *  35°. 
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Figure  12.  Yaw  vs.  Time.  With  No  Controller. 

UCK  =  18  Knots.  Rudder  Ordered  *  35°. 
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Figure  13.  Rudder  Response  vs.  Time. 

UCK  *  18  Knots.  Rudder  Ordered  *  35°. 
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Figure  14.  Deoth  vs.  Time.  With  No  Controller. 

UCK  »  12  Knots.  Rudder  Ordered  =  35°. 
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X  SCALE  =  40  Seconds  Per  Inch 


Y  SCALE  =  4  Degrees  Per  Inch. 


Figure  15.  Pitch  vs.  Time.  With  No  Controller. 

UCK  »  12  Knots.  Rudder  Ordered  =  35°. 
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Fiqure  17.  Speed  Change  vs.  Time.  With  No  Controller. 

UCK  =  12  Knots.  Rudder  Ordered  =  35°. 
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Figure  18.  Yaw  vs.  Time.  With  No  Controller. 

UCK  =  12  Knots.  Rudder  Ordered  =  35°. 
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Figure  19.  Rudder  Response  vs.  Time.  With  No  Controller 
UCK  ■  12  Knots.  Rudder  Ordered  *  35°. 
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Figure  20.  Depth  vs.  Time.  With  No  Controller. 

UCK  *  6  Knots.  Rudder  Ordered  =  35°. 
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Figure  22.  Roll  vs.  Time.  With  No  Controller. 

UCK  *  6  Knots.  Rudder  Ordered  *  35° 
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X  SCALE  =  40  Seconds  Per  Inch. 
Y  SCALE  =  40  Degrees  Per  Inch. 
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Figure  24.  Yaw  vs.  Time.  With  Wo  Controller. 

UCK  =  6  Knots.  Rudder  Ordered  *  35°. 
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IV.  AUTOMATIC  DEPTH  AND  PITCH  CONTROL 


A,  OPTIMAL  SOLUTION  TO  THE  LINEARIZED  SUBMARINE  EQUATIONS 
The  necessity  of  controlling  the  depth  and  pitch  of  a 
turning  submarine ,  as  was  stated  before,  stems  from  the  coupling 
between  the  states  such  that  any  changes  of  one  of  them  directly 
affects  the  other  one.  For  this  reason,  depth  and  pitch  control 
to  some  extent  must  be  accomplished  so  that  it  can  provide  a 
good  stage  for  the  roll  control.  There  are  many  ways  of  con¬ 
trolling  the  depth  and  pitch  which  depend  upon  the  way  of  using 
stern  and  fairwater  control  surfaces  combination.  The  scheme 


that  was  used  in  this  thesis  was  an  optimal  control  way  which 
was  originated  by  Drureys  in  Reference  7.  To  preclude  unneces¬ 
sary  repetition,  only  guidelines  of  the  method  and  the  differ¬ 
ences  of  this  study  are  to  be  discussed  here.  More  pronounced 
knowledge  can  be  found  at  Reference  7  and  8. 

It  can  be  seen  in  Appendix  A  that  the  equations  of  motion 
are  nonlinear.  To  be  able  to  use  this  set  of  equations,  they 
mu3t  be  linearized.  This  linearization  was  done  in  Reference  7, 
and  with  one  control  input  (sternplane)  they  can  be  written  in 
the  form  of:  _ 
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If  they  are  put  in  state  variable  form 
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The  problem  was  thought  as  a  linear  regulator  problem.  The 
general  scheme  of  linear  regulator  problem  is  shown  in  Figure 
26. 


The  cost  function  is 


tf 

7  J  [>*»  +  uTru1 


dt 


E  ■  X-r  and  r  is  the  reference  vector  which  is 


ordered  depth 

0 

r  * 

ordered  pitch 

a m 

0 

— 

49 


PLANT 


Figure  26.  Plant  And  Optimal  Feedback  Controller  For  Linear  Regulator  Problem 


As  a  result 


E  =  X  can  be  written. 
Following  the  Reference  7  and  8 


K  ■  -KA+KBR"1BTK-Q-ATK  and 
—  1  T 

U  =  -R  B  KE  can  be  written. 

If  the  linear  equations  are  augmented  and  put  in  the  state 
eguation  form  with  the  definition 
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In  this  study,  in  order  to  make  the  depth  rate  a  state  vari¬ 
able  the  following  trans Formation  was  used 

Depth  rate  *  w-u*  e  where  w  represent  the  comoonent  of 
u  in  the  z  direction,  Q  represents  the  nitch  angle.  In 


Reference  7  this  transformation  was  originally 

Depth  Rate  =  w  and  the  differences  between  these  two 
trams formations  showed  itself  in  the  calculation  of  K  values 
which  is  discussed  in  the  following  pages. 

Referring  to  the  cost  function 
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The  selection  of  the  weighting  factors  was  a  trial  and 

error  process.  The  series  of  weighting  factors  depending  upon 

the  relative  severity  of  the  influence  of  the  states  to  each 

other  was  tested  and  as  a  result 

E  «  Depth  Error  Weighting  =  10 

D  =*  Pitch  Error  Weighting  »  8000 

C  ®  Control  InDut  Weighting  =  100 

A  =  B  =  0  was  chosen. 

After  solution  of  ntn+l),  differential  eauations, 

'2 

K  values  associated  with  the  feedback  crain  was  found.  It  is 
observed  that  steadstate  values  of  K  are  constant  which  is  a 
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convenient  condition  for  ootimality.  Results  of  this  solution 
are  shown  in  Figures  27  through  36.  For  comparison  the  original 
A  matrix  of  Reference  7  with  the  transformation  Depth  rate  =  w 
was  also  used  to  calculate  K  values.  One  of  the  results  is 
shown  in  Figure  37.  It  is  observed  that  it  could  not  reach 
steady  state  value  in  the  reasonable  time  period. 

After  K  values  associated  with  the  feedback  gains  were 
found,  necessary  gains  were  found  via  the  optimal  low. 

1  T 

u  »  B  K  E 

The  result  was 

u  ■  DSAD  =  -0. 31623*ZOER-1.792*ZODOT  + 

36 . 069*PERR  =  102. 63*PlD0T. 

All  of  the  feedback  gains  were  found  at  15  knots  speed. 

It  is  obvious  that  these  gains  associated  with  the  hydro- 
dynamic  coefficient  used  are  function  of  the  speed.  To  make 
the  gains  compatible  with  the  soeed  range  of  the  ship,  they 
must  be  scaled  with  the  function  of  the  speed.  In  Reference  2, 
the  controller  gains  adjustment  as  a  function  of  the  speed  was 
discussed  deeply  and  found  that  gains  associated  with  the  depth 
rate  and  pitch  role  error  channel  was  inverselv  proportional  to 
the  instantaneous  speed  as  gains  associated  with  the  depth 
and  pitch  error  channel  remains  same.  After  the  justification 
of  these  results 

DSAD  *  -0.31623*Z0ER-1.792*(15)  .  *Z0D0T  + 

'UK 

36.069*PERR  +  102.63M15),  *PlD0T 

'UK 


3  3 


20.00  40.00  60.00  80.00  100.00  120.00  140 


20.00 


X  SCALE  =  20  Seconds  Per  Inch. 


Figure  31.  K22  vs.  Time.  UCK  =  15  Knots. 

D  »  8000,  A  -  0,  E  »  10,  B  =  0,  C  *  100. 
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Figure  36.  K44  vs.  Time.  UCK  *  15  Knots. 

D  a  8000,  A  a  0,  E  a  10,  B  *  0,  C  =  100. 
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X  SCALE  =*  20  Seconds  Per  Inch 
Y  SCALE  *  200  Per  Inch. 


Figure  37.  Kll  vs.  Time.  UCI<  *  15  Knots. 

With  Transformation  Of  De->th  Rate  »  W.  ' 
D  -  8000,  A  -  0,  E  *  10,  B  *  0,  C  -  100. 


DSAD  =  -0.31623n0ER-26.38.  *20D0T  + 

'UK 

36. 069*PERR  +  1539.45.  *P1D0T 

'UK 

where  UK  is  the  instantaneous  speed  (knots). 

B.  AUTOMATIC  DEPTH-PITCH  CONTROLLER  AND  SIMULATION  OP  THE 

RESULTS 

The  feedback  gains ,  which  are  the  ootimal  solution  of 
the  linearized  vertical  plane  equations  must  be  put  into  the 
general  controller  schema  in  which  the  submarine  dynamics  are 
represented  in  six  degrees  of  freedom  (i.e.,  non-linear 
equations  of  motion) .  Because  only  the  simulation  of  this  can 
justify  the  validity  of  the  solution.  The  complete  modified 
version  of  the  Drurey's  depth  and  pitch  controller  is  shown  in 
Fichre  38.  As  it  is  shown  the  sternplane  controller  (actuator) 
was  not  designed  as  a  part  of  the  controller  and  put  into  the 
controller  separately  with  the  plane  rate  7°/sec.  The  dynamics 
of  the  actuator  is  the  same  as  Drurey's  and  Stamps'  controller 
used.  The  same  actuator  dynamics  also  was  used  in  the  roll 
controller  which  is  to  be  discussed  in  the  following  section. 
The  actuator  block  diagram  is  shown  in  Figure  39.  To  preclude 
stability  which  occurs  as  a  result  of  the  excessive  uses  of 
the  sternplane  and  avoid  use  of  the  control  surface  with  big 
deflections  the  limiter  was  put  into  the  pitch  and  depth  error 
channels.  Pitch  error  limiter  was  10°  and  depth  error  limiter 
was  20  feet  and  found  by  a  trial  and  error  process. 

After  the  completion  of  the  controller  design,  various 
tests  were  run  to  see  the  effect  of  the  pitch  and  deoth  con¬ 
troller  on  depth,  pitch,  roll  and  sternplane  responses. 
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The  results  of  tests  at  6,  12,  18,  and  24  knots  to  a  constant 
35°  left  rudder  angle  are  shown  in  Figures  40  through  55. 

When  these  results  are  compared  with  the  responses  of  the  ship 
with  no  pitch  and  depth  controller  which  is  shown  in  Figures 
3  through  6  at  24  knots  and  in  Figures  8  through  25  at  18,  12, 
and  6  knots.  The  following  comparisons  can  be  drawn: 

1.  At  12  and  18  knots  controller  kept  the  ship's  depth 
and  pitch  stable  after  very  reasonable  transient  time. 

The  steady  state  pitch  and  depth  are  (5.2)°  -  (9.2)  feet 
at  12  knots  and  (4.7)°  -  (9)  feet  at  18  knots. 

2.  At  6  knots,  the  controller  failed  to  keep  the  pitch  and 
depth  stable.  The  ship  very  slowly  kept  on  losing  depth. 
The  reason  was  insufficient  speed  scaling  of  the  feedback 
gains  and  needs  for  another  control  surface  (Fairwater 
planes) .  Lack  of  fairwater  plane  control  surface  in  the 
depth-pitch  controller  unabled  it  to  keep  ship  stable  at 
low  speeds.  But  in  this  thesis  the  primary  concern  was 
to  minimize  the  snap  roll  which  occurs  at  high  speed. 

In  Reference  2  Stamps  collected  data  which  shows  the  oeak 
snap  roll  as  a  function  of  the  approach  speed  and  ordered 
rudder  angle.  This  data  is  repeated  at  the  following 
page  for  convenience  (Table  I).  It  is  seen  that  dangerous 
snap  roll  starts  occuring  at  approach  speeds  above  12 
knots.  For  this  reason  unability  of  the  controller  at 
low  speed  was  neglected.  Because,  a  control  scheme,  which 
at  low  speed  uses  Drury's  original  depth-pitch  controller 
(Fairwater  plane  is  part  of  controller)  and  at  high 


68 


speed  uses  this  design  (^airwater  plane  is  not  a  part 
of  depth-pitch  controller  but  is  a  part  of  the  roll  con¬ 
troller)  can  be  designed  and  switching  from  one  to  another 
can  be  achieved. 

3.  At  24  knots  the  controller  again  failed  to  keep  the  ship 
stable.  Depth  and  pitch  response  went  into  oscillation 
with  big  amplitudes.  Original  roll  response  of  the  ship 
which  is  shown  in  Figure  5  was  destroyed  in  the  sense  of 
decreasing  amplitudes  of  roll  oscillation.  It  did  not 
reach  stable  value  and  appeared  to  be  oscillating  in  the 
range  of  10°.  The  reason  was  the  following.  The  big  snap 
roll  which  was  around  37°,  initially  gave  very  big  dis¬ 
turbance  and  oscillation.  With  these  big  disturbances 
the  controller  which  uses  only  the  sternolane  as  a  con¬ 
trol  surface  was  unable  to  do  the  job.  It  was  thought 
that,  if  there  had  been  any  controller  which  could  have 
decreased  the  snap  roll  (Roll  controller)  it  would  have 
been  able  to  stabilize  the  ship  in  pitch  and  depth  as 
well  as  in  the  roll  response. 

Based  on  the  results  of  these  tests,  proceeding  with  the 
roll  controller  design  which  would  make  use  of  fairwater  plane 
as  a  control  surface  was  decided. 
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TABLE  I 


PEA*  SNAP  ROLL  ANGLE  (DEGREES)  VS.  APPROACH  SPEED  (KNOTS) 


AS  A  FUNCTION 

OF  ORDERED 

RUDDER 

ANGLE  (DEGREE) 

SPEED 

RUDDER 

(KNOTS)  4 

8 

12 

16 

20 

24 

1 

0.15 

0.  66 

1.45 

2.58 

4.03 

5.83 

2 

0.24 

0.99 

2.22 

3.95 

6.18 

8.88 

3 

0.31 

1.24 

2.80 

5.00 

7.84 

11.28 

4 

0.37 

1.45 

3.29 

5.87 

9.26 

13.27 

5 

0.41 

1.63 

3.71 

6.71 

10.32 

15.72 

6 

0.45 

1.80 

4.04 

7.39 

11.65 

17.90 

7 

0.49 

1.93 

4.46 

7.93 

13.11 

19.87 

8 

0.52 

2.  09 

4.81 

8.36 

14.47 

21.64 

9 

0.55 

2.22 

5.12 

9.10 

15.72 

23.23 

10 

0.58 

2.33 

5.38 

9.88 

16.88 

24.65 

11 

0.60 

2.42 

5.60 

10.62 

17.96 

25.94 

12 

0.63 

2.50 

5.78 

11.  31 

18.95 

27.13 

13 

0.65 

2.60 

5.94 

11.97 

19.86 

28.20 

14 

0.67 

2.72 

6.06 

12.58 

20.71 

24.17 

15 

0.68 

2.83 

6.33 

13.17 

21.49 

30.04 

16 

0.70 

2.92 

6.65 

13.72 

22.20 

30.83 

17 

0.72 

3.01 

6.96 

14.23 

22.87 

31.57 

18 

0.74 

3.09 

7.25 

14.72 

23.48 

32.22 

19 

0.76 

3.16 

7.54 

15.16 

24.04 

32.81 

20 

0.77 

3.22 

7.80 

15.59 

24.55 

33.  34 

21 

0.78 

3.27 

8.06 

15.98 

25.02 

33.82 

SPEED  (KNOTS) 
RUDDER 

4 

8 

12 

16 

20 

24 

22 

0.79 

3.32 

8.30 

16.34 

25.45 

34.25 

23 

0.80 

3.36 

8.52 

16.68 

25.84 

34.64 

24 

0.81 

3.39 

8.74 

16.99 

26.20 

34.99 

25 

0.82 

3.42 

8.95 

17.28 

26.51 

35.30 

26 

0.83 

3.44 

9.14 

17.55 

26.80 

35.57 

27 

0.84 

3.46 

9.32 

17.79 

27.06 

35.81 

28 

0.86 

3.47 

9.49 

18.01 

27.29 

36.02 

29 

0.87 

3.48 

9.65 

18.22 

27.50 

36.20 

30 

0.87 

3.49 

9.80 

18.40 

27.68 

36.36 

31 

0.88 

3.49 

9.94 

18.57 

27.  84 

36.50 

32 

0.89 

3.51 

10.06 

18.71 

27.97 

36.60 

33 

0.89 

3.58 

10.  18 

18.85 

28.09 

36.68 

34 

0.90 

3.64 

10.29 

18.96 

28.19 

36.75 

35 

0.90 

3.70 

10.39 

19.06 

28.27 

36.80 
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X  SCALE  »  40  Seconds  Per  Inch. 
Y  SCALE  =  5  Feet  Per  Inch. 


40.00  80.00  120.00  160.00  200.00 


Figure  40.  Depth  vs.  Time.  With  Deoth-Pitch  Controller 
UCK  *  6  Knots.  Rudder  Ordered  *  35°. 


X  SCALE  =40  Seconds  Per  Inch. 
Y  SCALE  =  2  Degrees  Per  Inch. 


Figure  41.  Pitch  vs.  Time.  With  Depth-Pitch  Controller. 
UCK  *  6  Knots.  Rudder  Ordered  =  35°. 
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Figure  42.  Roll  vs.  Time. 

UCK  »  6  Knots. 


With  Deoth-Pitch  Controller. 
Rudder  Ordered  *  35°. 


X  SCALE  =  40  Seconds  Per  Inch 


Y  SCALE  *  8  Degrees  Per  Inch. 


Figure  43.  Sternplane  Angle  vs.  Time.  T'Tith  Depth-Pitch 
Controller.  UCK  *  6  Knots  .  Rudder  Ordered 
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X  SCALE  =  40  Seconds  Per  Inch 


76 


o 


Figure  45.  Pitch  vs.  Time.  With  Depth-Pitch  Controller. 
UCK  »  16  Knots.  Rudder  Ordered  =>  35°. 


77 


,-12.00  -10.00  -8.00  -6.00  ->4.00  -2.00  0.00 


o 


X  SCALE  =*  40  Seconds  Per  Inch. 
Y  SCALE  ■  2  Decrees  Per  Inch. 
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Figure  46. 


Roll  vs.  Time.  With  Deoth-Pitch  Controller. 
UCK  *  12  Knots.  Rudder  Ordered  =  35°. 
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X  SCALE  =  40  Seconds  Per  Inch. 
Y  SCALE  *  2  Feet  Per  Inch. 
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Figure  48.  Deoth  vs.  Time. 

UCK  =  18  Knots. 


With  Deoth-Pitch  Controller 
Rudder  Ordered  »  35°. 
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X  SCALE  =  40  Seconds  Per  Inch 
Y  SCALE  =  8  Degrees  Per  Inch 


Figure  51.  Sternplane  Anqle  vs.  Time.  With  Depth-Pitch 
Controller.  UCK  »  18  Knots.  Rudder  Ordered 
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Figure  53.  Pitch  vs.  Time. 

UCK  -24  Knots. 


With  Depth-Pitch  Controller. 
Rudder  Ordered  *  35°. 


V.  AUTOMATIC  ROLL  CONTROLLER 


In  the  introduction  of  the  thesis,  the  possible  choices 
of  controlling  roll  was  discussed  briefly  and  stated  that 
fairwater  planes  were  to  be  used  as  control  fins  to  control 
the  roll.  The  reason  that  brought  uo  this  idea  lies  behind 
the  needs  of  more  control  surface  than  was  available  in  the 
present  design  criteria.  Roll  control  of  a  high  speed  sub¬ 
marine  has  continued  to  be  a  problem  over  the  years  because 
of  limited  control  surface.  Control  engineers  have  been  re¬ 
stricted  using  rudder,  stern  and  fairwater  plane.  All  of  the 
control  surface  was  not  meant  to  control  roll.  In  the  present 
design  criterias  stern  and  fairwater  planes  are  meant  to  con¬ 
trol  the  submarine  in  the  vertical  Diane  motion  (depth-pitch) . 
The  rudder  has  been  used  as  a  part  of  course  controlling  of 
the  ships.  Among  the  three  control  surfaces,  rudder  is  the 
only  one  that  has  direct  correlation  with  the  submarines  roll 
angles.  In  the  latest  study,  which  was  done  by  Stamps 
/Reference  27#  it  was  used  to  control  roll.  The  concept  of 
Stamps  roll  controller  was  based  on  the  idea  of  adjusting  the 
initial  rudder  angle  order  as  a  function  of  the  integral  of 
error  between  the  allowed  maximum  roll  and  actual  roll  ancles. 
The  initial  rudder  order  was  chosen  such  that  the  peak  roll 
expected  for  a  given  approach  speed  would  be  less  than  the 
maximum  allowed  roll.  The  integral  of  roll  error  was  then 
comDuted  and  scaled  to  represent  and  additive  term  applied  to 
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the  initial  rudder  order.  Stamps  design  didn't  give  any 
structural  changes  to  the  present  navy  submarines  design 
criteria.  In  a  sense  of  simplicity  it  was  perfect.  But  it 
prohibited  using  hard  and  excessive  rudder  order  which  is 
highly  desirable  at  some  submarines  required  maneuvering  in 
certain  tactical  areas.  Since  the  Stamos  design  slowed  down 
the  yaw  rate  and  caused  the  ship  to  change  its  course  very 
slowly,  alternative  design  ways  were  investigated  and  using 
fairwater  plane  in  the  differentially  deflected  mode  was 
chosen  as  a  possible  improvement. 

The  concept  of  roll  control  by  means  of  fairwater  planes 
is  based  on  the  idea  of  deflecting  the  fairwater  plane  dif¬ 
ferentially  such  that  it  can  give  a  roll  moment  in  opposite  direc¬ 
tion  to  the  instantaneous  roll  angle.  If  the  ship  has  a 
roll  angle  to  starboard  side,  the  Dianes  are  to  be  deflected 
to  give  a  roll  moment  to  portside.  The  positive  sense  of  this 
additional  roll  moment  created  by  the  differentially  deflected 
sailplanes  is  the  same  as  of  reference  5.  The  positive  sense 
of  sailplane  deflection  angle  adopted  in  this  thesis  (positive 
for  port  sailplane  deflected  leading  edge  up,  starboard  sail¬ 
plane  deflected  leading  edge-down)  is  in  agreement  with  that 
used  for  other  control  surfaces  in  Reference  5:  Positive 
when  the  surfaces  are  deflected  in  such  a  direction  as  to  in¬ 
crease  the  relevant  angle  of  the  submarine  about  its  mass 
center.  In  the  case  of  differentially  deflected  sailplanes, 
the  angle  is  the  roll  angle  and  this  is  defined  in  Reference  5 
to  be  positive  starboard  side  down. 
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Since,  in  the  present  design  criteria  of  navy  submarines, 
sailplanes  are  not  used  in  a  differentially  deflected  mode, 
standard  equations  of  motion  which  are  developed  by  NSRDC  in 
Reference  5,  do  not  have  correlative  terms  which  take  the 
additional  moment  term,  created  by  the  sailplanes  into  account. 
For  this  reason,  before  starting  to  design  the  roll  controller, 
this  additive  moment  must  be  estimated.  Once  the  counter 
moment  is  found  it  can  be  placed  in  the  righthand  side  of  the 
roll  moment  equation. 

A.  ESTIMATING  ROLL  MOMENT  DUE  TO  DIFFERENTIALLY  DEFLECTED 

SAILPLANES 

In  vertical  plane  motion,  the  normal  force  due  to  the 
fairwater  plants  (sailplanes)  is  given  by  the  equation  of 

fn  -  ^  L2u2-Z^  •  Sb 

Where  J  »  2 

L  =*  ship  length  in  feet 
U  =  forward  speed 

z6b  :  hydrodynamic  coefficient  associated  with  the 
sailplane  deflection  (in  conventional  mode  -  right  and  left 
side  moves  together  in  the  same  direction) . 

6b  ■  Deflection  of  the  sailplane  in  radians. 

After  setting J>»2  the  normal  force  equation  becomes 
fn  -  L2*U2*z£fcSb 

This  force  is  due  to  both  sides  of  the  sailplanes.  The 
force  due  to  one  side  of  the  plane  pair  is  half  of  the  total 
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1  2 

force  that  is  ^  L  1 


U 


Z6b&i? 


If  both  sides  of  the  planes 


are  deflected  by  the  same  amount  and  in  opposite  direction 
compared  to  each  other,  this  configuration  creates  a  moment 
associated  with  the  moment  arm  between  these  opposite  forces. 
This  is  shown  in  Figure  56.  Total  amount  of  the  moment  due  to 
sailplanes  can  be  written 


M 

M 


2*l/L2,U2‘Zg^  *Sb-  (moment  arm) 

L2* U2, Zgy  Sb.  (moment  arm) 
L3*U2*  Sb*  (Zg^  ‘moment  arm)yL 


if  this  momemt  term  is  placed  in  the  right  side  of  the 
equation  of  motion  about  the  body  axis  system  X-Axis  (roll 
axis) ,  the  complete  equation  takes  the  form  of 


Jxp  +  (Ix"IY)<Jr  = 


L5  /“Ko'  d  +  Kqr’qr  +  Kr'r  +  K_ 


>loj  plp!-7 

+  /~Kp'up  +  Kr'ur  +  K^,'v  +  Kwp  wo_7 

+  L3  /  *~K  *  *  u 2  +  Kv*uv  +  Kvlvl'v|(v2  +  w2)(1/2_7 


+  L3  K. 


+  L3-u2‘K5r‘6r 


3  2 

+  L  *U  • (Zr. ‘moment  arm) 


+  Sin  Cos 

In  this  last  form  Sb  is  replaced  by£>p  to  indicate  that  the 
fairwater  plane  is  to  be  used  in  differentially  deflected  mode 
(  ^^represents  deflection  anqle  of  fairwater  nlane  in  con¬ 
ventional  usage.  Both  right  and  left  sides  move  in  the  same 
direction). 
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Pigure  56.1.  Positive  Direction  Of  Differentially 
Deflected  Sailplanes. 

.2.  Negative  Direction  Of  Rolling  Moment. 


The  moment  arm  used  in  this  study  was  estimated  by  in¬ 
voking  the  data  from  Reference  6  that  belongs  to  the  un¬ 
classified  fictitious  submarine.  In  the  reasonable  prooortion 
to  the  length  of  the  ship,  the  moment  arm  was  chosen  as  7  ft. 
and,  according  the  data,  it  was  thought  acceptable. 

If  the  new  hydrodynamic  coefficient  is  defined  associated 
with  the  fairwater  planes  in  differentially  deflected  node 

z£f  -  (zSb  •moment  arm) 

°  L 

ZS?  »  (0.00558  x  7^251.75 

Z£*F  -  0.0015517. 

After  the  estimation  of  Z£*p,  as  a  hydrodynamic  co¬ 
efficient,  it  is  placed  in  the  modified  equation  discussed 
above.  After  these  modifications,  in  the  equations  of  motion 
all  of  the  term  which  is  a  function  of  the  5b  (sailplane  de¬ 
flection  in  the  conventional  mode)  was  set  to  zero  by  defining 
&>-  0.  As  result,  the  equation  of  motion  about  the  X-axis 
(roll  axis)  has  correlative  term  between  the  fairwater  plane 
and  roll  angle  of  the  ship  so  that  this  additive  moment  term 
can  be  used  to  give  feedback  to  smooth  the  snap  roll. 

B.  ROLL  CONTROLLER  DESIGN 

In  contrast  to  the  deoth  and  pitch  controller,  the  roll 
controller  was  designed  by  using  the  nonlinear  equations  of 
motion  and  linearizing  was  not  attempted.  In  the  previous 
section,  the  modification  of  the  roll  moment  equation  due  to 
the  imposition  of  moment  term  which  stems  from  the  fairwater 
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planes  deflection  was  discussed.  All  through  the  roll  con¬ 
troller  design  this  modified  version  of  the  equations  was 
used.  The  reason  that  led  to  using  six  nonlinear  equations 
was  such  that  linearizing  of  this  equations  in  six  dereee  of 
freedom  was  found  difficult  because  of  the  terms  coupled  with 
each  other  and  mainly  was  thought  that  from  the  output  of  the 
system  dynamics  the  roll  angle  and  the  roll  rate  (p)  could 
give  enough  feedback  information  to  accomplish  the  compensation. 

The  principles  for  designing  the  roll  controller  were 
simple.  If  the  ship  has  a  positive  roll  angle  (starboard) 
the  fairwater  plane  is  to  be  deflected  so  that  it  can  give 
negative  roll  moment,  i.e. ,  starboard  sailplane  is  to  be  de¬ 
flected  leading  edge  up  while  port  sailplane  is  deflected 
leading  edge  down  or  vice  versa.  As  long  as  the  ship  has  a 
roll  angle  this  would  cause  the  fairwater  plane  deflection  by 
way  of  the  feedback  channel  used.  In  the  design  it  was  assumed 
that,  the  fairwater  planes  actuater  was  capable  of  giving  e- 
qual  amounts  of  deflection  command  in  opposite  directions  to 
both  starboard  and  port  sailplanes.  Since  the  design  was  to 
invoke  extensive  simulation  study  (since  nonlinear  equations 
were  used)  in  determination  of  the  feedback  parameters,  the 
first  basic  controller  attempt  was  the  proportional  controller 
because  of  its  simplicity  to  design  and  implement.  The  pro¬ 
portional  controller  is  shown  in  Figure  57.1.  Referring  to 
the  figure,  if  the  submarine  has  a  roll  angle  this  would  give 
the  position  feedback  to  force  the  system  to  reach  zero  roll 
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Figure  57.1.  Position  Feedback  Controller. 

.2.  Position  And  Velocity  Feedback  Controller. 
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angle  position.  Determination  of  Kl  (proportional  constant) 
was  a  trial  and  error  process.  Even  though  position  feed¬ 
back  alone  gave  big  improvement  in  the  system  dynamics,  as  is 
to  be  discussed  in  the  following  pages,  it  failed  to  stabil¬ 
ize  the  system  in  some  operating  conditions  that  are  consider¬ 
ed  very  likely  to  be  faced.  This  leads  to  the  compensation 
of  the  system  by  velocity  feedback  and  other  modifications 
(Limiter  in  the  position  feedback  channel).  In  the  following 
pages  the  design  procedure  which  uses  extensive  computer  sim¬ 
ulation  is  discussed  from  the  simplest  case  of  proportional 
controller  to  the  last  modification  that  stabilized  the  system 
in  various  operating  conditions  and  over  a  wide  range  of  speeds 
1.  Proportional  Controller 

The  first  attempt  to  control  the  roll  was  to  design 
a  proportional  controller  as  is  shown  in  Figure  57.1.  Refer¬ 
ring  to  the  figure,  it  is  seen  that  ordered  fairwater  plane 
deflection  is  a  function  of  the  roll  angle  and  the  proportional 
constant  Kl.  Such  that 

ORDERED  FAIRWATER  DEFLECTION  *  DFOD  =  -Kl*ROLL 
Since  the  reference  signal  is  zero  the  system  always  looks  for 
zero  roll  angle.  Ordered  fairwater  deflection  as  an  input  to 
the  fairwater  planes  actuator  causes  deflected  planes  at  the 
output  due  to  the  actuator  dynamics.  Actuator  dynamics  were 
the  same  as  those  of  the  sternplanes  actuator  shown  in  Figure 
39.  With  the  counter  roll  moment  created  differentially  de¬ 
flected  fairwater  planes,  system  dynamics  try  to  decrease  the 
roll  angle. 


Determination  of  K1  was  a  trial  and  error  process. 

The  stability  range  was  1  through  4.  By  inspection  of  the 
results  of  computer  simulations,  K1  =  3  was  chosen  as  a  best 
choice  in  the  sense  of  snap  roll  and  steady  state  value  of 
the  roll  responses.  With  the  determined  Kl  =  3  value  the 
system  was  tested  at  24  knots  approach  speed  to  a  35°  rudder 
command  and  the  results  are  shown  in  Figures  58  through  62. 
Before  analyzing  the  results  it  is  instructive  to  indicate 
here  that  all  through  the  design  procedure  the  worst  condition 
that  could  happen  was  always  taken  into  account.  For  this 
reason  the  controller  was  designed  at  24  knots  base  speed  and 
to  a  35°  constant  rudder  angle  which  could  give  the  worst  snap 
roll  in  the  range  of  soeed  of  interest.  By  lookino  at  the 
Figures  58  through  62,  the  following  results  can  be  summarized: 

1.  As  was  predicted,  the  snap  roll  decreased  to  almost 

4°. 6  from  37°  and  roll  response  reached  a  steady  state 
value  of  3°. 5.  In  Figures  52,  53,  and  54  it  was  shown 
that  before  implementation  of  the  roll  controller  the 
system  was  unstable  to  an  identical  test.  Adding  the  roll 
controller  made  the  system  stable  at  high  speed.  In  depth 
and  pitch  responses  great  imorovement  has  been  made. 

Steady  state  values  of  the  depth  and  pitch  were  almost 
8  feet  and  4°. 5  which  were  acceptable,  while  this  deoth, 
pitch  and  roll  control  imorovement  has  been  obtained,  the 
stern  and  sailplane  was  used  moderately  such  that  they 
didn't  reach  any  saturation  (±  35°). 
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Figure  59.  Pitch  vs.  Time.  Kl  »  3.  UCK  =  24  Knots. 

Rudder  Ordered  *  35°.  Initial  Roll  Angle 
»  0°. 
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Figure  60. 


Roll  vs.  Time.  Kl  -  3.  UCK  =  24  Knots.  Rudder 
Ordered  =  35°.  Initial  Roll  Angle  =  0°. 
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X  SCALE  *  40  Seconds  Per  Inch 


Figure  61.  Sternplane  Angle  vs.  Time.  Kl  =■  3.  UCK  » 
24  Knots.  Rudder  Ordered  =  35°.  Initial 
Roll  Angle  *  0°. 
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Figure  64. 


Pitch  vs.  Time.  Kl  a  3.  UCK  =  24  Knots. 
Rudder  Ordered  .«  35°.  Initial  Roll  Angle 
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2.  Investigation  of  the  roll  response  in  Figure  60  lead 


to  further  study.  In  the  steady  state  value  of  the  roll, 

the  oscillation  with  almost  fixed  amplitude  and  fre- 

o 

quency  was  observed.  The  amplitude  was  around  0*1.  As 
a  first  guess,  the  oscillation  was  thought  to  be  a  stable 
limitcycle.  To  expose  the  problem  the  system  was  tested 
with  initial  roll  response  (5°  inboard  roll  angle)  at  24 
knots  to  a  35°  constant  rudder  angle.  Figures  63,  64, 
and  65  record  the  depth,  pitch,  and  roll  responses.  The 
system  became  unstable.  The  system  failed  to  compensate 
itself  when  it  was  commanded  to  turn  with  the  35°  rudder 
deflection  and  5°  initial  inboard  roll.  Any  submarine 
maneuvering  in  the  tactical  area  may  experience  a  turn 
with  initial  roll  angle.  Disturbances  from  the  heavy  sea 
state  might  cause  the  submarine  to  roll  eventhough  it  is 
in  a  straight  course.  If  the  submarine  is  commanded  to 
turn  at  this  moment,  it  was  shown  that  its  control  is  lost. 
That  is  why  the  situation  was  thought  unacceptable  and  com¬ 
pensation  of  the  system  was  attempted.  It  was  thought  that 
the  reason  for  the  failure  was  lack  of  enough  feedback  in¬ 
formation  and  compensation  of  the  system  with  velocity 
feedback  was  attempted. 

The  following  section  describes  the  modification  and  improve¬ 
ment  obtained. 

2.  Compensation  Of  The  System  With  Velocity  Feedback 
The  modified  controller  block  diagram  is  shown  in 
Figure  57.2.,  and  DFOD  *  -Kl*ROLL-K2*P. 
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It  was  felt  that  the  velocity  feedback  should  give  to  the 
system  better  dampino  and  better  stability  characteristics. 
Determination  of  the  K2  value  was  again  a  trial  and  error 
process  and  the  stability  ranqe  was  found  to  be  1  to  15  by 
computer  simulation.  By  inspection  of  the  results  of  the 
tests  with  different  K2  values,  K2  =  10  was  considered  the 
best  choice  in  the  sense  of  snap  roll  (max  roll)  and  steady 
state  value  of  the  roll  responses.  Since  the  reason  for  intro- 
dncincr  velocity  feedback  was  instability  in  the  presence 
of  an  initial  roll  angle,  tests  with  exaggerated  initial  roll 
(20°  inboard)  were  made.  It  was  Dreviously  denoted  that  the 
ship  has  almost  37°  snap  roll  at  24  knots  to  a  35°  rudder 
command.  To  force  the  ship  to  turn  with  a  hard  rudder  command 
when  it  already  has  a  very  serious  initial  roll  angle  was 
thought  a  good  example  of  the  capability  to  control  the  shin 
in  three  dimensions.  In  fact,  any  submarine  cruising  in  a 
straight  course  is  unlikely  to  face  this  much  roll  angle  from 
heavy  seastate.  In  Figure  66  the  result  of  this  test  is  shown. 
Before  the  implementation  of  the  velocity  feedback  the  system 
was  unstable  to  a  turn  command  with  initial  5°  roll  angle. 

And  now  it  is  stable  even  with  20°  initial  inboard  roll.  The 
result  was  such  that  steady  state  roll  angle  was  almost  4°. 5 
and  before  reaching  steady  value  or.  oscillation  with  decreasing 
amplitude  was  observed.  In  all  of  the  aforementioned  tests, 
the  ship  response  was  investigated  to  a  steady  35°  rudder 
command.  The  rudder  was  commanded  to  a  35°  deflection  at  the 
beginning  of  the  simulation  and  after  reachina  35°  it  was  not 
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X  SCALE  *  40  Seconds  Per  Inch. 
Y  SCALE  *  4  Degrees  Per  Inch. 


Figure  66.  Roll  vs.  Time.  K1  -  3,  K2  -  10.  UCK  »  24  Knots. 

Rudder  Ordered  =  35°.  Initial  Roll  Angle  »  -20°. 
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changed  to  any  other  command.  But  any  submarine  maneuvering 
in  tactical  areas  can  experience  successive  rudder  commands 
in  opposite  directions.  For  this  reason  the  test,  in  which 

Q 

the  rudder  deflection  was  commanded  to  35  in  the  time  inter¬ 
val  of  20  to  80  seconds  was  required.  Figures  67  through  71 
record  the  depth,  pitch,  roll,  sternolane  and  sailplane  de¬ 
flection.  In  this  test  the  initial  roll  angle  was  again  20° 
inboard.  The  result  was  unstable.  Inspection  of  these 
figures  reveals  a  very  important  reason  of  the  failure.  Curve 
number  2  in  these  figures  represents  the  rudder-  response.  But 
since  in  the  vertical  axis  the  automatic  scaling  was  used 
associated  with  the  output  responses  of  interest,  in  some 
figures  the  rudder  response  in  the  time  interval  of  2C  to  80 
seconds  was  seen  less  than  35°.  For  this  reason,  to  overcome 
misinterpretation,  curve  number  2  should  be  interpreted  as  the 
time  interval  where  rudder  deflection  was  35°.  By  inspection 
of  the  figures  67  and  68  its  seen  that  the  depth  and  pitch  re¬ 
sponse  are  almost  unchaged.  The  very  small  changes  are  due  to 
the  roll  oscilation  which  stems  from  the  initial  roll  angle. 

In  Figure  70  it  is  shown  that  in  the  time  interval  of  0  -  20 
seconds  the  sternplanes  oscillate  with  very  small  amplitudes 
to  compensate  these  depth  and  pitch  changes.  But  in  the  same 
time  interval  (0  -  20  seconds)  the  sailplanes  oscillate  with 
big  amplitudes  to  overcome  the  roll  oscillation  started  with 
the  initial  20°  inboard  roll  angle  (Figure  71).  Because  of 
the  roll  controller,  the  initial  roll  angle  has  more  effect 
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Figure  67.1.  Roll  vs.  Time.  Ki  *  3,  K2 

Knots.  Initial  Roll  Angle  =  -20 
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.2.  Rudder  Response  vs.  Time.  (Rudder  Ordered  *  35 
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Figure  68.1.  Pitch  vs.  Time.  Kl  »  3,  K2  *  10. 

24  Knots.  Initial  Roll  Ancle  =  -20' 

.2.  Rudder  Response  vs.  Time.  (Rudder  Ordered 
-  35°). 
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Figure  69.1.  Roll  vs.  Time.  Kl  *  3,  K2  »  10.  UCK  * 
24  Knots.  Initial  Roll  Anale  =  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered 
»  35°. ) 


12 


T 
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.2.  Rudder  Resoonse  vs.  Time  (Rudder  Ordered  = 
35°). 
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on  the  sailplanes  than  on  the  sternplanes.  In  the  tine  inter¬ 
val  of  0  -  20  seconds  of  Figure  69  the  roll  oscillation  started 
with  the  initial  roll  angle  is  seen.  Since  the  direct  effect 
of  the  roll  angle  on  the  sailplane  deflection,  this  initial 
roll  oscillation  causes  the  sailplane  to  oscillate  and  this 
oscillation  causes  initial  roll  oscillation  not  to  be  damned 
out.  As  this  is  going  on,  at  time  equal  to  20  seconds  the 
rudder  was  commanded  to  35°  full  deflection  and  this  gives  more 
unstability.  After  analyzing  this  and  previous  results  the 
following  conclusion  has  been  reached.  Since  the  roll  angle 
was  directly  fedback  via  Kl,  the  sudden  and  big  roll  response 
changes  cause  the  sailplane  to  oscillate  with  big  amplitudes 
and  this  leads  to  unstability.  To  overcome  this  difficulty  a 
limiter  was  placed  in  the  roll  error  feedback  channel  with  the 
magnitude  of  ±5  (The  same  reason  had  lead  to  the  placing 
depth  and  pitch  error  limiter  in  the  depth  and  pitch  controller 
design  as  was  discussed  in  Section  IV.B.).  After  this  last 
modification  the  system  was  tested  under  various  conditions 
which  are  to  be  discussed  below.  The  complete  controller  block 
diagram  with  the  depth  and  pitch  controller  is  shown  in  Figure 
72.  By  inspection,  the  simulation  results  of  the  last  version 
of  the  controller  can  be  sumarized  as  follows: 

1.  In  Figures  73  through  87  the  depth,  pitch,  roll,  stern- 
plane  and  sailplanes  responses  at  24,  18,  and  12  knots  to 
35°  rudder  command  in  the  time  interval  of  20  -  80  seconds) 
with  the  initial  20°  inboard  roll  angle  are  shown.  The 
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Figure  72.  Comolete  Depth-Pitch-Roll  Controller  (Final  Scheme) 


Figure  73.1.  Depth  vs.  Time.  Final  Result  With  Roll 
Error  Limiter.  Kl  *  3,  K2  *  10.  UCK  =* 

24  Knots.  Initial  Roll  Angle  -  -20°. 

.2.  Rudder  ResDonse  vs.  Time.  (Rudder  Ordered 
■  35°). 


117 


X  SCALE  =  40  Seconds  Per  Inch 


Limiter.  Kl  =•  3,  K2  »  10.  UCK  »  24  Knots. 
Initial  Roll  Angle  ■  -20°. 


.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  = 
35°). 
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Figure  75.1.  Roll  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  Kl  =  3,  K2  -  10.  UCK  =  24  Knots. 
Initial  Roll  £ngle  =  -20°. 

.2.  Rudder  Response  vs.  Time.  (Rudder  Ordered 
=  35°) 
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Figure  76.1.  Sternplane  Angle  vs.  Time.  Final  Result  With  Roll 
Error  Limiter.  Kl  =  3,  K2  =  10.  UCK  =  24  Knots. 
Initial  Roll  Angle  =  -20°. 

.2.  Rudder  Resoonse  vs.  Time  (Rudder  Ordered  =  35°). 
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Figure  77.1.  Sailplane  Angle  vs.  Time.  Final  Result  with  Roll 
Error  Limiter.  Kl  =  3,  K 2  =  10.  UCK  =  24  Knots. 
Initial  Roll  Angle  =  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  *  35°). 
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Figure  78.1.  Depth  vs.  Time.  Final  Result  With  Roll  Error  Limiter 
K1  a  3,  K2  *  10.  UCK  =  18  Knots.  Initial  Roll 
Angle  *  -20°. 

.2.  Rudder  Resoonse  vs.  Time  (Rudder  Ordered  =  35°) 
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Figure  79.1.  Pitch  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  K1  =  3,  K2  =  10.  UCK  =  18  Knots. 
Initial  Roll  Angle  =  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  =  35°). 
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Figure  80.1.  Roll  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  Kl  =  3,  K2  =  10.  UCK  =  18  Knots. 
Initial  Roll  Angle  =  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  =  35°). 
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Figure  81.1.  Sternplane  Angle  vs.  Time.  Final  Result  With 
Roll  Error  Limiter.  Kl  =  3,  K2  *  10.  UCK  = 

18  Knots.  Initial  Roll  Angle  =  -20°. 

.2.  Rudder  ResDonse  vs.  Time  (Rudder  Ordered  *  35°). 
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Error  Limiter.  Kl  =  3,  K2  =  10.  lTCK  =  18  Knots. 
Initial  Roll  Angle  =*  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  »  35°). 
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Figure  83.1.  Depth  vs.  Time  Final  Result  With  Poll  Error 
Limiter.  K1  -  >  K2  «  10.  UCK  ■  12  Knots. 

Initial  Roll  Angle  -  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  «  35° 
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Figure  85.1.  Poll  vs.  Time.  Final  Result  With  Poll  Error  Limiter 
Kl  =  3,  K2  »  10.  UCK  =  12  Knots.  Initial  Poll 
Angle  =  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  =»  35°). 
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Figure  86.1.  Sternplane  Angle  vs.  Time.  Final  Result  With 
Roll  Error  Limiter.  Kl  =  3,  K2  =  10.  UCK  * 

12  Knots.  Initial  Roll  Angle  =  -20°. 

.2.  Rudder  Resnonse  vs.  Time  (Rudder  Ordered  =  35°). 
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Figure  87.1.  Sailplane  vs.  Time.  Final  Result  With  ^oll  Error 
Limiter.  Kl  *  3,  K2  =  10.  UCK  *  12  Knots. 
Initial  Poll  Angle  =  -20°. 

.2.  Rudder  Response  vs.  Time  (Rudder  Ordered  =  35°). 


0.00 


Fiqure  88.  Depth  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  Kl  *  3,  K2  =  10.  UCK  =  24  Knots. 
Rudder  Ordered  =  35°.  Initial  Roll  Angle  =*  0°. 


132 


X  SCALE  ==.  40  Seconds  Per  Inch. 
Y  SCALE  =  .80  Degree  Per  Inch. 
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Figure  89.  Pitch  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  Kl  =  3,  K2  =  10.  UCK  =»  24  Knots. 
Rudder  Ordered  =  35°.  Initial  Roll  Ancle  *  0°. 
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X  SCALE  =  40  Seconds  Per  Inch. 


Figure  90.  Roll  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  K1  =  3,  K2  =  10.  UCK  =  24  Knots. 
Rudder  Ordered  =  35°.  Initial  Roll  Anale  =  0 
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Fiqure  91.  Sternplane  Angle  vs.  Time.  Final  Result  With 
Roll  Error  Limiter.  Kl  =  3,  K2  =  10.  UCK  = 

24  Knots.  Rudder  Ordered  =  35°.  Initial  Roll 
Angle  =  0°. 
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Figure  93.  Depth  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  Kl  =  3,  K2  =  10.  UCK  =  18  Knots. 
Rudder  Ordered  =  35°.  Initial  Roll  Ancle  =  0°. 
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Figure  96.  Sternplane  Angle  vs.  Time.  Final  Result  With 
Roll  Error  Limiter.  Kl  *  3,  K2  =  10.  UCK  * 

18  Knots.  Rudder  Ordered  ■  35°.  Initial  Roll 
Angle  ■  0°. 
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Figure  100.  Poll  vs.  Time.  Final  Result  With  Roll  Error 
Limiter.  K1  *  3,  K2  -  10.  UCK  *  12  Knots. 
Rudder  Ordered  =  35°.  Initial  Roll  Angle  =*  0° 
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X  SCALE  =  40  Seconds  Per  Inch 
Y  SCALE  =  4  Degrees  Per  Inch. 


Figure  101.  Sternplane  Angle  vs.  Time.  Final  Result  With 
Roll  Error  Limiter.  Kl  »  3,  K2  =  10.  UCK  * 

12  Knots.  Rudder  Ordered  =  35°.  Initial  Poll 
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same  test  without  roll  limiter  as  was  discussed  above 
gave  an  unstable  system  but  after  inserting  the  limiter 
the  result  was  stable  over  the  speed  range  of  interest. 

In  contrast  to  the  results  shown  in  Figures  67  through  71, 
there  is  no  oscillation  in  the  roll  and  sailplanes  re¬ 
sponses  in  the  time  interval  of  0  -  20  seconds  and  it 
was  observed  that  the  limiter  gave  very  good  damping  to 
the  system.  This  last  statement  can  be  justified  by 
inspection  of  Figure  75  such  that  at  time  equal  20 
(when  the  rudder  !•  commanded  to  full  deflection)  the 
roll  response  is  almost  zero  without  any  Drevious  oscil¬ 
lation.  As  the  stabilization  in  three  dimensions  was  being 
reached  the  stern  and  sailplanes  never  went  into  saturation. 

2.  Figures  88  through  102  record  the  depth,  pitch,  roll, 
sternplane  and  sailplanes  responses  at  24,  18,  and  12 
knots  to  35°  rudder  command  (the  rudder  is  commanded  to 
full  deflection  of  the  beginning  of  the  simulation)  with 
zero  initial  roll  angle.  It  is  observed  that  without 
using  excessive  stern  and  sailplane  deflection  great  sta¬ 
bility  in  the  three  dimension  has  been  reached.  As  was 
shown  before.  Figures  44  through  55  record  the  responses 
of  the  system  without  roll  controller  to  the  identical 
test.  Comparison  of  these  tests  is  made  in  the  Table  II. 

It  is  seen  that  at  high  speed  great  improvement  in  the 
sense  of  maximum  and  steady  state  value  has  been  obtained. 
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TABLE  II 


IMPROVEMENT  DUE  TO  ROLL  CONTROLLER 


WITHOUT  ROLL  CONTROLLER 

WITH  ROLL  CONTROLLER 

MAX 

STEADY  STATE 

MAX 

Depth  Changed 

8  ft. 

24 

Pitch 

UNSTABLE 

4°.  4 

Roll 

4° 

Depth  Changed 

10.1 

• 

00 

8.5 

8  ft. 

18 

Pitch 

9.5 

5° 

V£> 

• 

o 

in 

4°.  7 

Roll 

23° 

6° 

4°.  6 

3°. 2 

-  -  -  - 

Depth  Changed 

9.7  ft, 

9.5  ft. 

9.2 

9 

12 

Pitch 

7° 

5°.  7 

r- 

• 

to 

5.2 

Roll 

10°.  3 

2.5 

3°.  5 

2° 
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3.  Figures  103  through  105  record  the  roll  responses  of  the 
system  at  24,  18,  and  12  knots  to  35°  rudder  deflection 
(the  turn  was  commanded  at  the  beginning  of  the  simu¬ 
lation)  with  20°  inboard  roll  angle.  The  results  are 
stable.  Comparison  of  Figure  103  with  Figure  66  shows 
the  improvement  obtained  by  inserting  a  roll  limiter  such 
that  it  damped  out  the  oscillation. 

It  was  shown  that  utilizing  the  fairwater  planes  as 
the  control  surface  of  the  roll  controller  not  only  stabilized 
the  roll  response  but  also  gave  big  improvement  in  the  depth 
and  pitch  controller.  In  the  following  section  stability 
tests  of  the  system  to  disturbance  moments  is  discussed. 

C.  ROLL  CONTROLLER  STABILITY  TESTS 

Since  the  roll  controller  design  was  carried  out  usina 
nonlinear  equations  in  six  degrees  of  freedom,  the  stability 
analysis  of  the  system  was  not  practical.  To  investigate  the 
stability  computer  simulation  was  used  bv  applying  disturbance 
moments.  The  magnitude  of  the  moments  was 
DISTURBANCE  =  KS*(UCK)2*L3  where 
KS  »  4.0*0.0003*576/ (UCK) 2 

The  disturbance  moment  was  calculated  by  assigning  non¬ 
zero  value  of  KS  which  was  normally  zero  and  represents  the 
hydrodynamic  coefficient  which  gives  the  rolling  moment  when 
body  angle  and  control  surface  angles  are  zero.  The  magnitude 
of  the  disturbance  was  four  times  bigger  than  that  used  in  the 
stability  analysis  of  Stamps  controller.  At  24  knots  the 
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magnitude  of  the  disturbance  was  1.1028x10^  ft-lb.  such  that 
the  Stamps  controller  was  unable  to  stabilize  the  system  to 
the  disturbances  of  the  magnitude  in  the  order  of  10^.  Figures 
106  through  120  record  the  depth.  Ditch,  roll,  sternplane,  and 
sailplanes  responses  to  the  inboard  roll  moment  dis¬ 
turbance  mentioned  above.  The  distrubance  was  applied  at  time 
equal  20  as  a  step  moment.  Curve  number  2  represents  the  dis¬ 
turbances.  Inspection  of  the  results  shows  that  the  system  is 
stable  in  three  dimensions  to  the  disturbances  moment.  The 
worst  affect  of  the  disturbances  was  seen  at  low  speed  test 
which  was  normal.  But  even  in  this  case,  after  a  transient  re¬ 
sponse  the  system  was  stable. 
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Depth  vs.  Time.  Roll  Stability  Steo  Test. 

UCK  =  24  Knots.  Rudder  Ordered  »  35°. 
Disturbance  Moment  vs.  Time  (1. I028*107?t-lb. ) 
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igure  109.1.  Sternnlane  Angle  vs.  Time.  Roll  Stability  Steo 
Test.  UCK  *  24  Knots.  Rudder  Ordered  »  35°. 
.2.  Disturbance  Moment  vs.  Time  (1. 1028xl0^ft-lb. ) . 
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Figure  112.1.  Pitch  vs.  Time.  Roll  Stability  Steo  Test. 

UCK  =>  18  Knots.  Rudder  Ordered  =  35°. 

.2.  Disturbance  Moment  vs.  Time.  (1. 1028xl07ft-lb. ) 


Fiqure  113.1.  Roll  vs.  Time.  Roll  Stability  S 
UCK  =  18  Knots.  Rudder  Ordered 
.2.  Disturbance  Moment  vs.  Time.  (1 
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Figure  119.1.  Sternplane  Angle  vs.  Tine.  Poll  Stabilitv 
Step  Test.  UCK  =  6  Knots.  Rudder  Ordered 


.2.  Disturbance  Moment  vs.  Time.  (1.1028x10  ft- lb.) 
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VI.  DISCUSSION ,  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  WORK 

A.  DISCUSSION  AND  CONCLUSION 

The  characteristic  of  a  submarine  which  is  totally 
submerged  and  in  a  high  speed  turn  was  discussed  in  Chapter  II 
and  the  simulation  results  were  given  in  Chapter  III.  The 
unhydrodynamic  body  structure  when  it  is  submerged,  which  is 
mainly  due  to  the  appendages  (sail) ,  was  considered  the  main 
source  of  the  problem.  The  necessity  of  the  sail  that  provides 
rooms  for  periscope  and  other  vital  controls  that  are  located 
in  the  tower  does  not  permit  a  small  sail  structure.  On  this 
account,  a  control  scheme  must  be  used  to  compensate  the  prob¬ 
lems  that  stem  from  the  sail.  The  complexity  of  controlling  a 
submarine  in  a  hich  speed  turn,  comes  from  the  coupling  between 
states  and  three-dimension  dynamics.  In  Chapter  IV  an  automatic 
depth  and  pitch  controller  was  presented  briefly  with  the  modi¬ 
fication  of  the  original  design  represented  in  Reference  7. 

The  automatic  control  scheme  gave  the  advantage  of  taking  depth 
and  pitch  error  and  their  rate  into  account.  But,  since  only 
the  sternplane  was  used  as  a  control  surface,  the  results  made 
the  roll  controller  design  necessary  for  controlling  the  ship 
within  the  wide  range  of  speeds  and  under  various  disturbance 
conditions.  The  imposition  of  the  fairwater  olane  to  control 
roll,  which  was  the  main  goal  of  this  thesis,  gave  a  great 
amount  of  expected  improvement  to  the  system.  The  roll 


controller  scheme  was  carried  out  by  using  equations  of  motion 
(nonlinear)  in  which  the  vertical  plane  dynamics  were  stabil¬ 
ized  by  the  imposition  of  the  depth-pitch  controller,  and  roll 
state  and  roll  rate  gave  enough  information  to  provide  closed- 
Iood  control  for  the  ship.  The  control  scheme  was  a  direct 
roll  controller  compared  to  the  controller  designed  by  Stamps 
in  Reference  2.  Stamps’s  controller  was  an  indirect  roll  con¬ 
troller  and  in  his  design  the  rudder  was  used  as  a  function  of 
the  integral  of  roll  error  which  was  the  difference  between  the 
allowed  maximum  roll  and  the  actual  instantaneous  roll  angles. 
The  initial  rudder  ordered  was  to  be  such  that  the  peak  roll 
expected  for  a  given  approach  speed  would  be  less  than  the 
maximum  allowed  roll.  The  integral  of  roll  error  was  then  to 
be  computed  and  scaled  to  represent  an  additive  tern  applied 
to  the  initial  rudder  order  (9ro) .  The  instantaneous  rudder 
order  (DRCOM)  was  then  to  be  computed 
URCOM  -  Sro  *  Kr  (jZ*max  -0act>  at. 

The  technique  reported  in  this  thesis  has  the  benefits  of 
simplicity  and  costs  much  less  than  alterations  of  the  design 
criteria  of  Wavy  submarines.  But  since  it  prohibited  use  of 
a  hard  rudder  deflection  (35°)  in  the  early  phase  of  turn,  the 
time  required  to  pass  any  specified  yaw  angle  was  quite  long. 

For  this  reason,  the  roll  control,  in  a  direct  sense,  by  making 
use  of  fairwater  planes  as  control  fins,  was  investigated  in 
the  proposed  design.  As  a  result,  not  only  were  the  dangerous 
snap  rolls  decreased  to  magnitudes  of  5°,  but  also  with  the  hard 
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rudder  deflection  a  big  yaw  rate  was  achieved.  In  Table  III, 
the  performance  of  the  submarine  with  each  of  these  two  roll 
controllers  is  compared.  In  the  Starnes*  design,  maximum  roll 
angle  is  predetermined  and  the  controller  is  designed  not  to 
exceed  this  value.  The  higher  allowed  maximum  roll  angle  gives 
higher  yaw  rates.  In  the  proDosed  desion,  maximum  roll  angle 
is  the  result  of  the  system's  dynamics. 

TABLE  III 

COMPARISON  OF  THE  PROPOSED  AND  STAMPS ' DESIGN 


TYPE  OF 

MAX 

APPROACH 

TIME  REQUIRED  IN 

SECONDS  TO 

CONTROLLER 

RQLL 

SPEED 

PASS  SPECIFIED  YAW  ANGLE 

■ 

45° 

90° 

180° 

STAMPS'  DESIGN 

Mi 

24 

52 

91 

18 

39 

67 

■H 

12 

27 

47 

BUM 

DIFFERENTIAL 

5°. 61 

24 

11.  5 

20 

41 

FAIRWATER 

4°.  69 

18 

14 

25 

54 

PLANE 

3°.  4 

12 

19 

36 

80 

It  is  seen  that  the  use  of  differential  fairwater  planes 
provides  a  faster  response  with  reduced  maximum  roll  angle  at 
all  soeeds. 

In  the  proposed  design,  only  one  fairwater  nlane  actu¬ 
ator  was  used  and  it  was  assumed  that  it  gives  deflection  in  a 
differential  mode.  Switching  criteria  from  the  conventional 
usage  of  the  fairwater  nlane  (both  starboard  and  port  sides 
move  simultaneously  in  the  same  direction)  to  the  differential 
mode  was  not  discussed. 
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This  research  has  shown  that  imposition  of  the  fairwater 
plane  as  a  part  of  the  roll  controller  improved  the  control  of 
a  turning  high  speed  submarine  in  three  dimensions.  Dangerous 
snap  roll  was  decreased  to  around  5°  and  with  only  sternplane 
and  the  stabilized  roll  response  high  improvement  in  the  depth 
and  pitch  characteristic  was  achieved. 

B.  RECOMMENDATIONS  FOR  FURTHER  WORK 

The  following  research  would  be  worthwhile  in  future 
studies : 

1.  Automatic  depth  and  pitch  controller  by  using  only 
sternplane  in  the  sense  of  an  optimal  feedback  controller 
for  linear  tracking  should  be  designed  and  the  combina¬ 
tion  of  this  design  with  the  present  roll  controller 
should  be  investigated.  This  study  would  give  good  in¬ 
sight  into  the  problem  which  comes  from  the  depth  changing 
when  the  submarine  is  in  a  high  speed  turn. 

2.  Switching  criteria  from  the  roll  controller  of  Reference 
2  to  the  proposed  roll  controller  should  be  established 
and  these  two  designs  should  combine  together.  This  would 
give  better  depth  and  pitch  control  to  the  ship  when  it 
has  low  speed  where  a  great  amount  of  roll  control  is 
not  needed. 

3.  With  the  coordination  of  NSRDC  a  more  accurate  hydrodynamic 
coefficient  associated  with  the  differentially  deflected 
sailplane  should  be  obtained. 


4.  Addition  of  an  integral  of  error  to  the  roll  control 
should  be  studied  to  see  if  the  steady  state  roll  angle 
can  be  reduced. 
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APPENDIX  A  • 

EQUATIONS  OP  MOTION 


The  following  set  of  equations  are  referred  to  a  body 
fixed  systcn  of  axes  which  arc  coincident  with  the  prln 
clpal  axes  of  Inertia  of  the  body.  The  origin  of  this 
axis-system  Is  located  at  the  assumed  center  of  mass  of 
ehe  bodv 

Equation  of  Motion  Along  the  Body  Axis  System  x-Axis 
m(i  -  vr  ♦  wq).  1 1*  [x  ^  •  q*  t  Xrr  ’  r1  +  X  rp  •  rp] 

* 

*  Oi  ■«♦*,,,  •»«■] 

♦f*‘  [x».  •»*] 

*  f ‘  V  [X  .rir  ’  »r  ‘  < *  X  *  X  ,Mb  '  «b'] 

*  -  ‘I"' 


♦  X  «.W  '  «*■  -  » 


sin  0 


Equation  of  Motion  Along  eh*  Bod/  Axis  System  y-Axta 


m(v  -  wp  ♦  ur)  *2  l*  £y  *  '  r  ♦  Y  ^  '  p  J 

♦I1*  [%q  ' p<*  *  yp!p!  *pW] 

♦*?**  [Y  v  '  *  *  Y  wp  '  WP  +  Yvjr|,|7|l<vl  +  w*>*  I  M] 
♦  |la  [Yr  ur  >  Y  |rj  jr  '  u|  r|  6r  ♦  Y  p  *up] 

*iX*  Yra'  '  <n*  *  l>ttp 


*fX'  TY*  '  tt*+  Yv  '  uv+  Yv|v|  '  v!(v2  ♦ 

w*>*|] 

tf*VY 

♦f*VY4rn*  ‘ 

♦f*‘Yvn.  '<«•-  »)uy 

Yv|v|n*  '  («'  -  D  v|(v*  ♦  w*)*  | 

*■§  i*  Y  »  wv# 

2  wv 

X*  (F„) _  V*  ♦  w*  (-w)  Sin.ut 

•  7  ▼  *  Q 

ff  Multiplied  by 

u 

U 

*  EW|  tin  4>  coa  0 

for  large  angle*  of 
attack  near  -90® 

Equation  of  Motion  Along  the  Body  Axis  System  t-Axls 


m{w  -  uq  *  vp)  a  |  <‘  Z  ^  ^ 

‘  f1*  [*„  •'*»♦*„  ••»•] 

Note  1 

♦  -£ A3  [Z  w  '  "  f  2  y,  '  vr#  +  2  yp  ,vP+*Zvp  '  VP»] 

♦  liJ  [z  ,  '  u*  +  z  |,|6.  '“Iq|«»*zw|q|'^|(v,  +  wl)*|  W] 

♦f4*  2  qn'  <«*  •  «  «q 

♦  f1*  [z  •  ’ttl+Z  W  'Uw  +  Zw|w|  'WI(V*  +W*)il] 

♦  2  t*  fZ|  w|‘  UM  +  2WW'  |w  (v*  t  w1)'!  ♦  Zyy'  v*»  J 


♦!iV  [Z  6s  ’4'  +  26b  ’4b] 


*  lX*  [2  wo*  '  {n'  -  l)  UW  +  Z  w|w|B,  «  (n'  .  1)  w|(v*  +  wl)i|J 


♦flV2W  (»*  -  1)  8# 


♦  -f  i*  (F ,)  v*  +  w1  y  Sin  wt 

a  »  va  q 

♦  EWj  cos  9  co.  <t>  ,  Multiplied  by  * 

u 

.  IT 

for  large  angles  of 
attack  near  -90* 

’  Note  I 

when  not  multiplied  by 
u  add  to  Z  • 

V  *P 
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Equation  o(  Motion  About  the  Body  Axis  System  x-Axle 


lj>  *  a,  -  ly)  qf  »  I  *•  [k  J  •  P  *  K  qr  •  qr  ♦  K  .  •  f  ♦  K 


p|p|  '  i 


[K,  '“P*Kr 

'ur  +  K;  '^+Kwp  '  wp 

-f*3 

[K.  '«‘.Ky 

'  uv  +  K  vjvj  *  v|(v*  ♦  w*)‘ 

**• 

K  1  vw 

yw 

f** 

“‘K  ftr  ’  4r 

♦  B*  B  sin  $  coe  8 


177 


Equation  of  Motion  About  the  Sod/  Axle  System  y-Axle 


Note  1 

Iy4  *  (1*  '  V  rp  *  iL%  [M  q  *  «  +  M  rr  ’  tU  +  M  rp  '  rP  *  &M  rp  '  rP*] 

♦  {l*  [Mq  '  “«I  ♦  M  |q|6a  '  u|qj«*  ♦  M  jw|q  '  |(v*f  w*^Jq] 

+  {t*  [Mw  '**Mvr  ’vr#+Mvp  ,Vp#] 

♦  li4Mqn'  ,(n*  '  *>“* 

+  -|ia  [m  „  '  u1  +  M  ^  •  uw  +  M  wj  wj  '  w|(v* +w,)^|J 

♦f*3  [M|w|’  u^wl  *  Mww’  lw  (v*  +  w*)*l  +  Mvv’  ] 

+  f**  u*  [m5<  •#.♦»«*  -Sb] 

+  f  /3  M  wn*  ’  tn*  "  11  uw 
+  IZ3  Mw|w|n'  '  ln'  "  »"l<v*  ♦*»>*! 


*Pa 

♦  Be  B  ain  0 

•  EW  j  co*  0  cos  $ 


#  Multiply  by  rj  tor 
large  ancles  or  attack 
near  -9Cr 

Note  I 

when  not  multiplied 
by  «  add  to  Mrj>* 
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f.'qua  t Ion  of  Motion  About  tha  Body  Axla  Systea  x-Axla 


J.rMJy  *Jx)pq  -|i‘  [n.  •  *♦  Npq  'p,+  'p*] 

♦  1 14  [n  r  •  ur  ♦  N  jr|4f  •  u|r|«r  ♦  N  jvjr  '  |(yl  ♦  w,)^|r] 

*it*  [Np  ‘up+n*  '  wp] 

+  -£x4  Nrn,  *  {*’  -  1)  ur 
« 

♦  |xa  [„.  *  «*♦  Nv  '  UV4-  NV|VJ  «  v|(v*  +  w1)*  |] 

♦  f.’  «*N4,  •«, 

♦|«*  “*Nlra'  '»«, 

♦  I  xa  N  vn,  »  (n’  -  1)  UY 

♦  li3  Ny|v|n'  ’  {n’  ■  »)vHv«+w*>f| 

♦  •£  ta  N  *  wVf 

2  wv 

♦  cos  9  »ln<J>  I  Multiply  by  j  for  large 

angles  of  attack  near  -90* 


AUXILIARY  EQUATIONS 


t 


<J)  ■  p  4  i  »ind 

•  •  .  . 

0  m  (q  -  (icos9  Sirup  )/  COI<P 

♦  ■  (r4  0sin4>)/cos9cos4> 

ico  a  ucos0cos$4v(  *in<t>  tin  6  cos  #  .  cos  <{)  sin$) 
4  w  (  sin<$>  sin  $  4  cos  <£>  sin  0  cos 

y«  a  u  cos  0  sin  i>  4  v  (  cos<(>  cos  $  4  sin$  sin  0  sin  0) 
4  w  (cos<j)  sin  9  si nj>  -  sin<£  cos  0) 


*0  * 

-  a 

sin  9  4  v 

cos  0  sin$  4  w  cos  0 

COS  4> 

• 

U  « 

(U* 

4  vJ  4  w 

(Fx)p 

S 

[  4  a,’  n’  4  n’*  ’ 

1  when 

kj  <  n' 

9 

f  X*«* 

[  V  4  V  ♦  bs’  " 

|  when 

k,  <  n'  <  k* 

9 

f  A*u*| 

[  •«,'  +  c,'  n'  4  c,'n'*  * 

|  when 

kj  <  n'  <  k. 

9 

[  d,'  4  V  n1  4  d3*  n'»  ] 

|  when 

n'  <  kg 

«  •  • 
•l*  %2* 

J  h. 

•Is  ®2»  b3 

•  i  . 
«1#  e2»  «3 
s  t  • 
*1#  *2» 


Sets  of  non-dl-.enstoi.nl  coefficients  used  In  ths  pro¬ 
pulsion  equation  above.  Th-  set  vr.icl'.  will  be  tn  effect 
st  any  tlxe  diirir.,-  a  slrvila^ci  rnneuver  will  depend  on  ths 
vsluo  of  n*  and  the  numbers 


NOMENCLATURE 


All  symbols  used  la  the  equations  of  notion  and  In  tho 
auxiliary  equations  and  relationships  which  appear  In  this 
report  ars  defined  below.  Any  dimensions  Involved  will  be 
consistent  w*ch  the  foot-pound-second  system  of  units.  All 
angles  are  In  degrees.  The  Foreran  variables  corresponding 
to  these  symbols  are  shown  In  Appendtx  g  . 


a 


SYMBOL 


S 


SL 

v 


u 


DEFINITION 

A  dot  over  any  symbol  signi¬ 
fies  differentiation  with  res¬ 
pect  to  time. 

Buoyancy  force  which  Is  posi¬ 
tive  upwards. 

Maes  of  the  submarine  Including 
the  water  In  the  free  flooding 

spaces. 

Overall  length  of  the  submarine 

Linear  velocity  of  origin  of 
body  axes  relative  to  an  earth- 
fixed  axis  system. 

Component  of  U  along  the  body 
x-axis. 

Component  of  U  along  the  body 
y-axts. 

Component  of  U  along  the  body 
t-axls. 
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Conn  and  sp««d:  A  steady  value 
of  u  for  a  given  propeller  rpm 
whence,  and  control  surface 
ingles  aro  aero.  Sign  changes 
with  propollsr  reversals 

Longitudinal  axis  of  the  body 
fixed  coordinate  axis  system. 

Transverse  axis  of  the  body 
fixed  coordinate  axle  system. 

Vertical  axis  of  the  body  fixed 
coordinate  axis  system. 

Distance  along  the  x0  axis  of  sa 
earth-fixed  axis  system. 

Distance  along  the  yQ  axle  of  an 
earth-fixed  axle  system. 

Diet*146*  along  the  »’0  axle  of  an 
earth-fixed  axle  system. 

Component  of  angular  velocity 
about  the  body  fixed  x-axls. 

Component  of  angular  velocity 
about  the  body  fixed  y-axls. 

Component  of  angular  velocity 
about  the  body  fixed  s-axls. 

• 

The  a4coordtnate  of  the  center 
of  buoyence  (CB)  of  the  subma- 


1 


- 

o< 

Angle  of  attack. 

13 

Annie  of  drift. 

• 

• 

OA 

c r 

V 

Dh 

Deflection  of  bowplane  (  or 
eallplane  ) 

Sr  , 

Dr 

Deflection  of  rudder. 

Ss, 

D* 

Deflection  of  scernplans. 

•* 

The  ratio  «c/u. 

ft 

m 

6 

Pitch  ancle. 

• 

$ 

Taw  ancle. 

* 

t 

*oll  ancle. 

m  * 

c 

Maas  density  of  sea  water. 

o 

tft 

WelRht  of  water  blown  from  a 
particular  ballast  tank  ldent- 
Ifled  by  the  Integer  assigned 
to  the  Index  1. 

6> 

Angular  velocity. 

e 

Time. 

*ti 

location  along  the  body  x~axta 

of  the  center  of  mass  of  th# 

ballast  tank  when  this  tank  is 
filled  with  sea  water. 
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* 


(r  )_  Propulsion  (ores  (so«  auxt- 

*  P 

llary  equations  and  relatlon- 
ships). 

I  Moment  of  Inertia  of  a  sub¬ 

marine  about  the  x-axis. 


Moment  of  Inertia  of  a  sub¬ 
marine  about  the  y-axls. 


*1 


Moment  of  Inertia  of  a  sub¬ 
marine  about  the  s-axls. 


KP-  V-  KP I  p|#  •*<»/. 
K*'.  ■'up'.  Ko'. 


K 


/ 

V 


*  ir  * 

"  Hr 


Non-dlmenslonal  constants  each 
of  which  Is  assigned  to  a  parti¬ 
cular  force  term  In  the  equation 
of  motion  about  the  body  x-axls. 


M. M  *,  M  AM  M  *  M.  i-V 
q  '  rr  *  rp  tp  q 


'Mq*  V*  Mvr  *  *  Mvp’’  K*’* 


V’  Mv|vr  M|w|’*  Mvw’*  Mvv** 
MAb'*  **vn *  *  *  Mwjw|n'  '  M4sn* 


Non-df mens  Iona  1  constants  eaeh 
of  which  is  assigned  to  a  parti¬ 
cular  force  term  In  the  equation 
of  motion  about  the  body  y-axls. 
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4M&.V  - 


V*  V*  V-  V*  VUc*-  *1  vl  r*  * 

V«  V*  Swp'*  ** rn *  *  *  V*  V* 


»vM*-  S«r’’  ^drn’ ' '  X  vo  * #  *  Mv|  v|  «• 


H  • 
w» 


Non-dlnens tonal  constants  tach 
of  which  Is  assigned  to  a  parti* 
eular  fore*  torn  In  tha  equation 
of  aotloa  about  tha  body  z-axts. 


i  f  r  1 
*  *cr  • 

V* 

V  *  V-  xw. 

t  l  y  • 

UU  •  *VV  • 

T  * 
WV  • 

X4  r4  r*  *  X4a4a' 

X  •  x 

4  bd  b  *  vvn 

»  f  t  Y  • 

•  wwa*  *  e«6an*  • 

*4  rd  rn '  * 

Non-dimensional  constants  aach 
of  which  la  assigned  to  a  parti* 
eular  force  term  In  the  equation 
of  notion  along  tha  body  x-axls» 


Tv|r|’*  V’  Y|r|4r’’  Yp '  •  Yrn’'1 
V*  V-  Yv|v|'.  Y4  r'  •  Y4rn' ' ’ 


Yvn  * '  *  Y 


Vl  V 


!«*’•  v-  <y 


vs 


Non-dinenstonai  constants  each 
of  which  Is  assigned  to  a  parti¬ 
cular  force  term  In  tha  equatlo* 
of  notion  along  the  body  y*axls 


AZvp’’  V*  z!«il«s*  zv|,r* 


1  1 
qn 

*n*.  V* 

Zw|w|’‘ 

Vi’* 

.  • 
wv  • 

Xvv’’  Z4.' 

'•  Z4b  *  * 

2wn ' ' * 

wlw| 

n”  Z4sn” 

<Vv. 

Non-dimensional  constants  each 
of  which  Is  assigned  to  a  parti¬ 
cular  force  tern  in  the  equatin' 
of  notion  along  tha  body  z*axls 
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